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We consider a two dimensional Turing like system with two diffusing species which 
interact with each other. Considering the species to be charged, we include the 
effect of an electric field along a given direction which can lead to a drift induced 
instability found by A.B.Rovinsky and M.Menzingeryj. This allows one to study 
the competition between diffusion and drift as was done numerically by Riaz et al. 
We show here that an analytic formula can be found on the basis of a linear stability 
analysis that incorporates all the effects that are known for the system and also 
allows for some detailed predictions. 

Pattern formation in two reacting and diffusing species was first studied by Turing 
Turing argued that if the diffusion coefficients of the two species are widely different, then 
if one of the species is auto catalytic with the other inhibiting its growth, then the steady 
homogeneous state will be unstable to a patterned steady state. The instability could also set 
in as a temporal pattern in a spatially homogeneous state under certain conditions. Turing 



patterns have been a very important aspect of the study o 



nonlinear systems 



mmm. 



^, Q, llll that pattern formation 



Decades later it was found by Rovinsky and Menzinger 
could occur even if the two diffusion coefficients were nearly equal provided there was an 
external electric field and the diffusion species were charged giving rise to a gradient coupling. 
The pattern formed in this case would be travelling waves as opposed to the standing waves 
of the Turing pattern. Recently there has been extension of the Turing work in some 
unexpected directions [l2,[l^. In a recent study, Riaz et alQ, 1^ have shown numerically 
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that a Turing pattern for charged species could be altered by an applied electric field [l6|. 

In this work, we revisit the pattern formation problem with an external electric field to 
arrive at a single analytic result that has the property of capturing all the possible conditions 
for the instabilities in the system. We work with a two dimensional set up as in the work 
of Riaz et al. We take the system to be unbounded in the y direction and be bounded by 
rigid plates at x = ±L. The boundary conditions are that the concentration of the species 
vanishes at the boundary and so does the current normal to the plates which is proportional 
to the x-derivative of the concentrations. The electric field is taken to be in the x direction 
which leads to a drift in that direction. The existence of the plates in the x-direction is 
vital to keep the system bounded. The plates also play the very important role of fixing the 
wave number in the x-direction. In the absence of the constraint there will be an overall 
selection mechanism for the wave numbers (k = ki + but the individual components 
are not uniquely determined. What we will see below is that ki is fixed by the boundary 
condition and thus once is known both A;iand k2 will be determined. In the work of Riaz 



et al 14| the fixing of ki and k2 is undertaken numerically. Here the analytic fixing of ki 



allows us to find a general expression for the thresholds of the different instabilities. 



The general reaction-diffusion prob 
modelled by the evalutaion equation [14 1 



em for two species A{x,y,t) and B{x,y,t) can be 



8 A f)A 

^ ^ DV'A + z,ED— + f(A.B) (1) 
§.V^B,.yji,,(A,B) (2) 

In the above, D is the diffusion coefficient for the species A in units where the diffusion co- 
efficient for the species B is unity. The external electric field in the x-direction in denoted by 
E and zi and Z2 are the charges associated with the inhibitor and the activitor respectively. 
The operator is two dimensional (pattern on a substrate) and the function f{A, B) and 
g{A, B) describe the growth and interaction of the species A and B. The electric field terms 
come from an expression for the current together with the relevant Einstein's relation. In 



the Gierer-Meinhardt model 



13, 
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f{AB) = ^-A + a (3) 
g{A,B) = f^{A'-B) (4) 

The growth rate of A due to interaction with the substrate is a and the natural decay rate 

□ 

for B is /i. In the Lengyel-Epstein model[19| 

A B 

f{A,B) = <yKB-^^) (5) 

g[A.B) = a-B--^^^ (6) 

where a, h and a are constants. The homogeneous steady state is A = Aq and B = Bq 
such that f{AQ,Bo) = qIAqjBq) = 0. The hnear stabihty analysis around A = Aq and 
B = Bq leads to 

^ = DV'SA + ^2ED^^ + anM + a.^SB (7) 
^ = V^^i? + z^E^^ + a^^SA + a,,5B (8) 

where an = (li)^^,^^, «i2 = (1^)^0,50' ^^i = (I|)a„Bo ^22 = {%) a,,b. 

We consider a geometry which is confined by plates at x = ±L and is unbounded in the 

?y-direction. The solution will be periodic in the y-direction and if we take the wavenumber 

in this direction to be ^2? then we can write 

(M, 6B) = {Ai{x),Bi{x)) exp {ik2y) exp (Xt) (9) 

where A is the eigenvalue determining the temporal growth. Then Ai{x) and Bi{x) satisfy 
the differential equation 

[A + Dkl - D— - Z2ED— - an]^i = a^^Bi (10) 

[A + - ^ - ziED— - 022]^! = 021^1 (11) 
Eliminating Bi one can write 



D ^+'^^+[iDa22 + an)-2Dkl-{l + D)X 



- X'^ -kl{au + Da22) + Dkl] = (12) 
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where 

= ED{zi + Z2) 
P2 = E[ziaii + Z2Da22 — {zi + Z2)Dkl] 
^oT" = On + 022 - (1 + D)kl 

A = 011O22 - 012O21 (13) 

At this point the general procedure should be clear. We need to solve the homogeneous 
4th order equation above. This will involve four arbitrary constants which have to be fixed 
by boundary conditions. Since the system is homogeneous the four conditions will lead 
to four homogeneous linear algebraic equations and for consistency the determinant has to 
vanish. The resulting equation fixes A in terms of L, k2, E and other system parameters. 
The requirement ReX > for instability allows us to discuss the different situation that can 
occur. The above procedure in general and in principle cumbersome. We illustrate this in 
the simpler situation of £" = and the Turing limit i.e. D << 1. The lesson that we learn 
here will be put to good use for the more complicated case. 

Accordingly, we set = in Eq.(12) and obtain 

+[A' - A + Dkl - kl{au + Da22) + A] A = (14) 
The general solution can be written as 

4 

Ai = ^ Q exp {jniix) (15) 

where j = \/— l.Then from eq.(14) we have 

Dm'^ - [Da22 + an - ^Dkl - A(l + D)]m'^ 
+ [A2 - '^A + Dk^ - kj{an + 022^) + A] = (16) 

For D << 1 (The Turing case) the two roots are approximately (we keep Dk^ since k2 is 
not known a-priori) 

mi - (17) 



D 

ml c - 7 1--^'^ (18) 



A^ - aA - k^aii + Dk^ + A 
an - A - 2DH 



5 



where a — an + 022 — ^l- 

For even solutions we can write 

Ai — ci cos(mix) + C2 cos(?ri2x) (19) 

Imposing the boundary conditions Ai = and no flux condition = on x = ±L one 
can write 

nil tan(miL) = m2 tan(m2l/) (20) 
In the hmit D << 1, mi — > 00 and this require m2-L — > | leads to the condition 

-[TrA - {kl + ^)]A - {kl + ^)an 
+ A + Dkt + 2Dkl— = Q (21) 

where TrA — an + 022- 

We note that + ^ enters as a combination which we call k'^ and to 0(D). we can add 
to eq.(21) a term D{^Y without committing any significant error since L is very large as 
well. In that case eq.(21) becomes 

A^-(rrA-A;2)A-anA;2 + L>A;^ + A = (22) 

This reproduces the Turing condition to the leading order in D, since we find the condition 
for instability is 

A-^<0 (23) 
and the characteristic wave number kmin is given by 



kLn = ^ (24) 



The lesson that we learn from the above exercise is that the operator ^ can be effectively 
replaced by and using eq.(12), we determine A from 

A^ - X[TrA-{l + D)k^ + —E{zi + Z2D)] 

ITT 

+ Dk" - k\an + a22D) + A - —EDz2{k^ - 022) 

1 

- i^Ezi{Dk^ - an) - j^ZiZ2E^D = (25) 
The above equation can be written as 

A^ - A(ai + ia2) + A + i/?2 = (26) 



6 



where 

ai = TrA - (1 + D)k'^ 

/3i = Dk^ - k\an + ^022) + A - —E^ziZ2D 

/?2 = ^^[-^l^ll + -^2^022 - (-21 + Z2)De] (27) 

The real part of the eigenvalue A in eq.(25) will be negative (the condition that the homo- 
geneous state will be stable) provided 

f3i < ai{aiPi + a2p2) (28) 
using «!, a2, Pi, P2 given in eq.(27) leads after straightforward algebra to the central result 

{zi - Dz2f[aiia22 - k\au + ^022) + Dk^]—E^ 

> [TrA - (1 + D)e]^[k\au + ^022) - A - Dk^] (29) 

This is the single formula that contains all the possibilites of the pattern formation in the 
presence of the electric field. 

We will now discuss the various possibilities. 

First consider the Turing problem E — 0. Now, the right hand side of eq.(29) has to be 
negative for stability which means 

A-k'^{au + Da22) + Dk^ >0 (30) 

for stabihty. This is the central criterion for stabihty for all A > in the presence of 
diffusion. If the sign is reversed in eq.(30), we get the Turing instability and for A > 0, this 
occurs for a band of wavenumber where the above expression is negative. The minimum of 
the expression is obtained for k'^ — ""20°^^ value at the minimum is A — ^"""'^^"^^^ 

and hence the instability criterion is A < ^"^^^^"^^'^ , an inequality which is easy to satisfy 
for D «1 

We now consider the opposite limit i.e. there is no diffusion and only drift. In this case 
eq.(29) acquires the form 

{zi - Dz2faiia22^E'' > (TrAfi-A) (31) 
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Since we want to start with an initially stable state i.e. TrA < and A > 0, we have 
< and eq.(31) becomes 

ttE \TrA\ f -A Y 

< I ' J I — (32) 



■kEq _ 


TrA 






2L 




-DZ2 




\a11a22J 



2L \zi - Dz2\ \a11a22. 
This clearly shows that an instability will set in if E > Eq, where 

(33) 

We see immediately that for the instability to set in one must have a differential mobility i.e. 
Zl 7^ Dz2. This result in accordance with A.B.Rovinsky and M.Menzinger 9|. In the case 
of D ~ 1, i.e. the two diffusivites are nearly equal( a situation very different from Turing), 
we get for instability 

7]-2 

{zi - Z2Y[aiia22 - k^{aii + 022) + ^1^-^^ 

> {TrA - 2A;2)2[A;2(an + 022) - A - A;^] (34) 

We treat the situation which for i? = is stable so far as the reaction goes and is also 
stable when diffusion is included. This implies A > 0, TrA < , A — (an + 022)^;^ + A;^ > 0. 
The right hand side of eq.(34) is now negative and for the inequality to hold, we need 
011022 — A;^(oii + 022) + With 011O22 < and TrA < 0, we can satisfy this inequality 

in the range < k < ko where 



2kl = TrA + v/(rrA)2 - 4011O22 (35) 

and for k < ko, the critical value of E which will trigger an instability will be given by 

TT^ {TrA~2eY e{au + a22)-A-k^ 

> —JZ .^2..^ T?in. , „ N , /.4 (36) 



4L2 [zi - Z2y 011022 - A;2(aii + 022) + k^ 

According to the above relation all wavenumbers greater than k^ are always stable. In 
writing down the above condition, we see the advantage of the exact expression of eq.(29). 
The order of magnitude estimation in 9|] does not yield the above answer. 

For any D, we note that if diffusion destabilizes a stable reactive system, then A;^(oii + 
-D022) — A — Dk^ > and it follows that 011022 — k\au + ^022) + Dk^ — O11O22 — A + 
[A — A;^(oii + -D022) + Dk^] < since 011O22 < 0. Hence eq.(29) can never be satisfied. No 
amount of electric field can stabilize the system. Finally, if the diffusive system is stable 
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FIG. 1: Plot of critical electric field {E) vs wave number (k) for an =0.899, 022 =-0.91, ai2 =1, 
021 =-0.899, zi=l and Z2 =2 

i.e. k^{aii + 1)022) — A — Dk'^ < 0, then a critical E{Ec) will destabilize the state provided 
«ii«22 — k'^iO'ii + Da22) + Dk'^ < which will happen \i k < k'^ given by 

, '2 _ Qii + -DQ22 + \/ (qii + Da22Y - 4:Daua22 , 
1^0 - ^ (>^7j 

For small enough D, k'^^ iaii+Da22) 

^ ^ [TrA-il + D)kr ( -A \ 

In summary we have studied a reaction-diffusion system in the presence of a constant 
electric field along a particular direction. We have found a single analytic expression which 
contains all possible information about the stability and the instability of the system for 
different ranges of the diffusion coeffecient. The primary result that emerges is that there is 
an uppar limit on the wave number of the instability and that for each wave number below 
that there is a critical electric field that can excite that particular wave number. 
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